Crumpling transition and flat phase of polymerized phantom membranes 
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Polymerized phantom membranes are revisited using a nonperturbative renormalization group 
approach. This allows one to investigate both the crumpling transition and the low-temperature, 
flat, phase in any internal dimension D and embedding dimension d, and to determine the lower 
critical dimension. The crumpling phase transition for physical membranes is found to be of second 
order within our approximation. A weak first-order behavior, as observed in recent Monte Carlo 
simulations, is however not excluded. 
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Membranes form a particularly rich and exciting do- 
main of statistical physics in which the interplay between 
two-dimensional geometry and thermal fluctuations has 
led to a lot of unexpected behaviors going from flat to 
tubular and glassy phases (see [1, 2, 3, 4] for reviews). 
Roughly speaking, membranes fall into two groups [4]: 
fluid membranes, in which the building monomers are 
free to diffuse. The connectivity is thus not fixed and the 
membrane displays a vanishing shear modulus. In con- 
trast, in polymerized membranes the monomers are tied 
together through a potential which leads to a fixed con- 
nectivity and to elastic forces. While fluid membranes arc 
always crumpled, polymerized membranes, due to their 
nontrivial elastic properties, exhibit a phase transition 
between a crumpled phase at high temperature and a 
flat phase at low temperature with orientational order 
between the normals of the membrane [4, 5, G, 7]. Amaz- 
ingly, due to the existence of long-range forces mediated 
by phonons, the correlation functions in the flat phase 
display a nontrivial infrared scaling behavior [8, 9, 10]. 
Accordingly, the lower critical dimension above which an 
order can develop appears to be smaller than 2 [10], in 
apparent violation of the Mermin- Wagner theorem. 

Let us consider the general case of D-dimensional 
non self-avoiding (phantom) membranes embedded in a 
d-dimensional space. Early e-expansion [7] performed 
at one-loop order on the Landau-Ginzburg- Wilson-type 
model relevant to study the crumpling transition of poly- 
merized membranes has led to predict that just below the 
upper critical dimension D ~ A, the crumpling transition 
is of second order for d > dcr = 219 while it is of first or- 
der for d < dcr ■ This leaves however open the question of 
the nature of the transition in the physical (D = 2, d = 3) 
situation, the case e = 2 being clearly out of reach of such 
a one-loop order computation. On the numerical side for- 
mer Monte Carlo (MC) studies (see [11, 12] for reviews) 
predict a second-order behavior while more recent simu- 
lations [13, 14] rather favor first-order behaviors. There 
is however no definite conclusion and no explanation for 
these versatile results. 

In parallel to the investigation of the crumpling tran- 



sition, an effective elastic field theory has been used to 
probe the flat, low-temperature, phase of membranes 
[4, -5, 8, 10]. An e-expansion has been performed [8], 
also at one-loop order, below the upper critical dimen- 
sion D = 4 showing that this flat phase is controlled by a 
nontrivial flxcd point (FP). However, again, this low or- 
der computation performed in the vicinity of D = 4 has 
been of no use to accurately determine the properties 
of genuine 2D membranes such as the critical exponents 
and the lower critical dimension Dic{d) above which the 
flat phase can exist. 

Signiflcant progress has been realized with the use of 
large-d expansion [10, 15], and variant of, such as self- 
consistent screening approximations (SCSA) [16] that 
have allowed to evaluate the exponent 77 both at the 
crumpling transition and in the fiat phase as well as the 
dimension Dic{d). However, the very nature of the ap- 
proach, requiring large values of d, makes doubtful the 
quantitative predictions extrapolated at small d and even 
impossible the determination of the line dcr{D), separat- 
ing the first-order from the second-order regions. 

A flaw of the previous approaches to polymerized mem- 
branes is that, due to their perturbative character, they 
arc unable to treat all aspects of the physics of mem- 
branes including crumpling transition, fiat phase, and 
lower critical dimension and thus to get a global picture 
of the renormalization group (RG) phase diagram. In 
this article, we propose an approach of polymerized mem- 
branes based on a nonperturbative RG method [17] which 
has been applied successfully in both particle and con- 
densed matter physics (see [18, 19, 20] for reviews). With 
this method, that we adapt to the treatment of extended 
objects, we are able to describe within the same formal- 
ism, i.e., using a unique effective action and a unique set 
of RG equations, both the crumpling transition and the 
flat phase of membranes. Concerning the crumpling tran- 
sition, we reproduce the results obtained within the e- 
expansion approach and at leading order within the large- 
d approaches. Moreover, we determine the line dcr{D) 
everywhere between D = A and D = 2. Our estimates of 
dcr{D = 2) ~ 2 lead to predict a second-order phase tran- 
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sition in the physical case but do not completely exclude 
a weak first-order behavior. Our investigation of the flat 
phase also allows to recover all previous perturbative re- 
sults including e and l/d-expansions. Moreover we get, 
for all values of d, a determination of the lower critical 
dimension Dic{d) above which the crumpling transition 
and flat phase fixed points are shown to coexist. 

Our approach is based on the concept of effective av- 
erage action [17] (see [18, 19, 20] for reviews), rfe[r], 
where r = r(x) is a d-dimensional external vector that 
describes the membrane in the embedding space while x 
is a set of internal ZJ-dimensional coordinates which la- 
bels a point within the membrane. The quantity rfe[r], 
k being a running scale going from a lattice scale fc = A 
to the infrared scale k ~ 0, has the physical meaning 
of a coarse grained free energy where only fluctuations 
with momenta q > k have been integrated out. Thus, at 
the lattice scale A, F^^a identifles with the continuum 
limit of some lattice Hamiltonian while at long distance, 
i.e., at A: = 0, it identifies with the standard free energy 
r. The fc-depcndence, RG flow, of Ffe is provided by an 
exact evolution equation [17], 
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where t = Ink /A. The trace has to be understood as a 
i'-dimensional momentum integral as well as a summa- 
tion over internal indices. In Eq.(l), Rk{q) is an effective 
infrared cut-off function which suppresses the propaga- 
tion of modes with momenta q < k and makes that Tk 
encodes only modes with momenta q > k. A conve- 
nient cut-ofltis provided by Rk{q) = Z{k^ - q^)9{k'^ - q^), 
where Z is a field renormalization - see below. It gen- 
eralizes to translational-invariant action density a cut-off 
[21] which has been largely used since it leads to com- 
pact expressions. Finally, and of utmost importance, 

(2) 

note that the term F^ in Eq.(l) is, in principle, the 
exact, i.e., full field-dependent, inverse propagator, the 
second derivative of F^ with respect to the field r, taken 
in a generic, nonvanishing field configuration. This is the 
fact at the very origin of the nonpcrturbative character 
of the method. 

Let us now make precise the form of F^. It must 
be invariant under the group of Euclidean displacements 
which includes translations and rotations. This imposes 
to Ffe to be a functional of daV = dr/dxa, a = 1 . . . D, 
the order parameter, and of scalars in both the embed- 
ding and membrane spaces. An exact treatment of Eq.(l) 
would imply F^ to enclose all powers and derivatives of 
these Euclidean invariants. This goal is however unrealis- 
tic and one has to truncate Ffe. We choose here an ansatz 
that allows both to make easily contact with previous - 
perturbative ~ approaches and to realize our program. It 



is given by: 

Ffc (r) = J '^"^ \ {dadarf + u[dar.di3r - C'^Sap)'^ 
+ v{d^Y.d^v-DCf (2) 

where Z, u, v and are the running couplings which 
parametrize the model, with the indices a and /3 run- 
ning over 1...D. This is, up to a redefinition of the cou- 
plings, the action used in [7] to investigate the crumpling 
transition. Let us recall the physics encoded in Eq.(2) 
at the mean-field level with u > and u + vD > 0. 
For = 0, the minimum of F^ is given by a configu- 
ration where daV vanishes which characterizes a crum- 
pled phase. For > this minimum is given by 
a configuration r(x) = C^a=i^a ^a, where the {e^} 
are D orthonormal vectors, which corresponds to a D- 
dimcnsional flat phase. Action (2) thus describes a tran- 
sition between a high-temperature, crumpling, phase and 
a low-temperature, flat, phase. The excitation spectrum 
in the ordered phase is provided by d — D out-of-plane, 
capillary, waves and D in-plane, phonon, modes. A cru- 
cial aspect of our approach is that, since we establish 
nonpcrturbative RG equations for the couplings enter- 
ing in Eq.(2), and in particular for the coupling ^, we are 
able to tackle both the crumpling transition, typically as- 
sociated to a vanishing ^, and the flat phase flxed point 
(FLFP) which is reached by letting ( run to inflnity. 

Technically, the flow equations for the couplings Z,u,v 
and ^ are obtained using their definitions in terms of 
functional derivatives of the effective action (see [19, 20, 
22] for details) and applying RG Eq.(l). In terms of 
dimcnsionlcss quantities, these equations write as follows: 
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D{D + 2y V ; on 
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+ {d-D) [u^ + 2{D + 2)uv + D{D + 2)v^)l^+' 
+ ((3D + 2)u^ + {D^ + D- 2){iuv + Dv^))lg+^ 
+ {9u^ + 6{D + 4)uv + {D^ + 6D + 12)v^)l^+'^} 

where Ad = 2-°~^n-^^'^ /T{D/2). The flow of Z, that 
provides the function rit = —dlnZ/dt giving the critical 
exponent ry at a FP, is too long to be displayed here (see 
[22]). In Eqs.(3) /fj,^ is a shortcut for : 
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where Pi{q) = Z + Rk{q) + ,i = 0,1,2 and 

d/dt only acts on Rk- These so-called "threshold func- 
tions" (see [19, 20]) control the relative role of the differ- 
ent modes, phonons and capillary waves, within the RG 
flow. In Eq.(4), the mass mo = is associated to the 
d — D transversal, capillary, modes while m\ = AC,'^u and 
m\ = 8^^ {u + v) are masses associated to the D phonons 
modes that split up into D — 1 modes with mass mi and 
one mode with mass m2. 

- The crumpling transition - Let us first consider the 
crumpling transition. To recover the RG equations de- 
rived pcrturbativcly in [7] one expands Eq.(3) in powers 
of both e = A — D and the couplings u and v that are 
of order e at any putative nontrivial FP. This also corre- 
sponds to an expansion in powers of the phonon masses, 
which are small at the crumpling transition FP. Using 
the fact that the threshold functions entering in the flow 
of u and v have a universal, cut-off independent, limit at 
vanishing masses in Z? = 4 given by l^^^^ ~ 1, one obtains: 
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Up to a change in variable {v v ~ u/A) these are 
the equations derived in [7]. We recall that, at suffi- 
ciently high values of d, i.e., d > dcr = 219, just be- 
low D = A, the sets of Eqs.(3) and (5) admit a sta- 
ble (in the u and v directions) FP associated to the 
crumpling transition, called crumpling transition fixed 
point (CTFP). Still at c? > dcr, there exists another FP, 
close to the CTFP, which is unstable and that, when 
the dimension d is lowered to dcr, annihilates with the 
CFTP, defining the curve dcr{D). A large-d analysis 
of Eqs.(3) can be also easily done. The leading contri- 
butions come from the capillary modes which enter in 
Eq.(3) through the terms proportional to d — D. With 
our cut-off function I^qq = A/D and — 
that the coordinates of the CTFP are given by (^r — 
16Ad/D{D^ - 4), Ucr (16 - D^)D{2 + D)/(256dA£,) 
and Vcr = -(16 - D^)D/{256dAD). The correspond- 
ing critical exponents arc: i/ ~ 1/{D — 2) + 0{l/d) and 
r] = 0{l/d) in agreement with [ti] and [10]. 

To tackle with the physics below D = 4 we have nu- 
merically solved the FP equations between D = A and 
D = 2, a dimension in which the effects of truncation 
start to be important. The right part of Fig.l summa- 
rizes our results: one finds a smooth curve dcr(D) which 
starts at dcr = 219 in D = A and reaches dcr — 2 in 
D = 2 leading to predict a second-order phase transition 
for physical membranes. In this last case one finds, at the 
CTFP, a thermal exponent v = 0.52 and rj = 0.627 which 
compares well with the results provided by the large-d ex- 
pansion T] = 2/3 [(i, 10] and MC results r] = 0.71(5) [23] 
but less with the Monte Carlo Renormalization Group 
r] = 0.85(15) [24] and the SCSA rj = 0.535 [1(.]. At our 



level of approximation, our results display a weak de- 
pendence with respect to the cut-off function Rk (q) that 
induces an error on the curve dcr{D). Using another 
cut-off, Rk{q) = Zq'^/{exp{q'^/k'^) — 1), we have evalu- 
ated the error bar on dcr{D ~ 2), which is typically of 
order 5dcr ^ 1- This means that one cannot exclude 
dcr[D = 2) to be close to, or even slightly above, d = 3 
so that the crumpling transition for genuine membranes 
would be predicted to be of weak first-order in agreement 
with recent MC results [13, 14]. This point will be further 
analyzed in the near future [22]. 
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FIG. 1: On the right part, the curve dcr{D) which separates 
the region with a CTFP and without a CTFP. On the left 
part, the lower critical dimension Dic{d). 

- The Rat phase - The equations relevant to study 
the flat phase are easily obtained in our formalism by 
considering the regime C » 1 in the RG flow Eqs.(3), 
which corresponds to a regime where the phonon masses 
are very large and thus to a regime dominated by the 
fluctuations of the capillary waves, as expected in the 
deep flat phase. Setting d = d — D one gets: 



dtu^{D-A + 2r]t) u + 
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{D^ + 6D^ + 8D'2){u + v) + 128{D^ - l)u{u + 2v)Ad 
with Ad = AoiS + D - rjt) and, for a = 



dta = {D — 2 + rjtja — 



16 d{Q + D - T]t)a'^ Ad 
D{D^ + 8D + 12) 



(7) 



an equation which generalizes, to any value of D and d, 
the one obtained in the limit of large elastic constants, 
D = 2 and large-d, in [G]. Note that the function rjt in 
Eqs.(6) and (7) determines, at a FP, the exponent rj of the 
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capillary waves. The analog exponent, r]u, for the phonon 
modes, is obtained by the usual Ward identity [10]: r]u ~ 
4 — D — 277 that follows from rotational invariance. 

The set of Eq.(6), when expanded in powers of e = 
4 — _D, degenerates into those derived perturbatively in 
[8]. Accordingly, Eqs.(6-7), admit three nontrivial FPs, 
among which one, the FLFP, is stable with respect to 
all directions including a down to a dimension Die, the 
lower critical dimension. In the limit of large codi- 
mension d the coordinates of the FLFP are a/ = 
iCj 00), Uf = (16 - D^)D{2 + D)/{256dAD) and 
Vf = —(16 — D^)D / {256dAD) , with these two last quan- 
tities being identical to those of the CTFP. At the FLFP 
one finds 77 = 0{l/d), in agreement with previous large-d 
approach [10]. Moreover Eq.(7) indicates that, at large 
d, the FLFP is stable down to Dic{d 00) = 2. in 
agreement with [10] which predicts: Dic{d 00) = 
2 - 2/d + 0(1/^2). Note also that the RG flow on u 
and V indicates that for d = and at any nontrivial FP 
one has the exact result ?] ~ {D — A)/2 [10]. For physi- 
cal membranes one finds j] — 0.849 which compares well 
with the SCSA 77 = 0.821 [10] and numerical simulations 
T? = 0.750(5) [23] and t] = 0.81(3) [25] but less with the 
large-d result 77 = 2/3 [10]. 

Finally Eqs.(6-7) also allow a determination of Dic{d) 
for all values of d. To do this we use the equality [10]: 
r]{Dic,Uf,Vf) = 2 — Die at the FLFP which defines Die 
as the dimension at which phonons and capillary waves 
scale identically [10]. In fact, the RG Eq.(7) provides 
another interpretation of Die- Indeed, above Die, Eq.(7) 
possesses a solution with a ^ which corresponds to the 
CTFP. Indeed, just above Die, one has D — 2 + ?]er «C 1 
so that the nontrivial solution of Eq.(7) obeys aer «C 1 
and thus (^er ^ 1 which is precisely the regime of validity 
of this equation. Thus, just above Die, Eq.(7) well de- 
scribes both the CTFP and the FLFP which coexist and 
Die corresponds to the dimension at which the CTFP 
collapses to the FLFP which becomes unstable. Using 
the relation r]{Die, uf,Vf) = 2 — Die one obtains: 
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which, once inversed, provides the expression of Die{d). 
The corresponding curve is displayed on the left part of 
Fig.l. In particular one has, for genuine membranes: 
Die{d = 3) ~ 1.33. This result displays a remarkable 
stability with respect to a change in the cut-off function. 
With Rk{q) = Zg^/(exp((jVfc4) - 1) one finds Die{d = 
3) ~ 1.30. Our results compare well with the large-d, 
Die{d = 3) = 4/3 [III], and SCSA, Die{d = 3) = 1.5 [IG]. 

In summary, we have investigated the crumpling tran- 
sition and flat phase of Z)-dimensional polymerized mem- 
branes embedded in a d-dimensional space within a non- 
pcrturbative RG approach. We have determined the 



whole line der{D) that separates the second- and first- 
order regions, the lower critical dimension Die(d) and the 
critical exponents. More sophisticated ansatz should be 
used to systematically increase the accuracy of our results 
although implying a heavy algebra. Finally, our approach 
can be applied to many other situations in which the per- 
turbative approaches lead to an unsatisfying quantitative 
or even qualitative description such as in self-avoiding, 
anisotropic or disordered membranes [1, 2, 26]. 

We thank F. David, B. Dclamotte, T.C. Lubensky, L. 
Radzihovsky and J. Vidal for helpful discussions. 



[2: 

[3: 

[4] 
[5] 

[6: 

[7 

[s: 

[9 

[10: 
[11 

[12 
[13 

[14 
[15 



(8) [16; 



[17 

[is: 

[19 

[20 

[21 
[22' 
[23^ 
[24; 

[25; 

[26; 



* Electronic address: kownacki@u-cergy.fr 
t Electronic address: mouhanna@lptmc.jussieu.fr 
[1] D. R. Nelson, T. Piran, and S. Weinberg, eds., Proceed- 
ings of the Fifth Jerusalem Winter School for Theoretical 
Physics (World Scientific, Singapore, 2004), 2nd ed. 
See the contribution of L. Radzihovsky in [1]. 
M. J. Bowick and A. Travesset, Phys. Rep. 344, 255 
(2001). 

See the contributions of D. R. Nelson in [1]. 

D. R. Nelson and L. Peliti, J. Phys. (Paris) 48, 1085 
(1987). 

F. David and E. Guitter, Europhys. Lett. 5, 709 (1988). 
M. Paczuski, M. Kardar, and D. R. Nelson, Phys. Rev. 
Lett. 60, 2638 (1988). 

J. A. Aronovitz and T. C. Lubensky, Phys. Rev. Lett. 
60, 2634 (1988). 

E. Guitter, F. David, S. Leibler, and L. Peliti, J. Phys. 
(Paris) 50, 1787 (1989). 

J. A. Aronovitz, L. Golubovic, and T. C. Lubensky, J. 
Phys. (Paris) 50, 609 (1989). 
See the contribution of Y. Kantor in [1] . 
See the contribution of G. Gompper and D. M. Kroll in 
[1]- 

J. -Ph. Kownacki and H. T. Diep, Phys. Rev. E 66, 
066105 (2002). 

H. Koibuchi et. al., Phys. Rev. E 69, 066139 (2004). 
M. Paczuski and M. Kardar, Phys. Rev. A 39, 6086 
(1989). 

P. Le Doussal and L. Radzihovsky, Phys. Rev. lett. 69, 
1209 (1992). 

C. Wetterich, Phys. Lett. B 301, 90 (1993). 

C. Bagnuls and C. Bervillier, Phys. Rep. 348, 91 (2001). 
J. Berges, N. Tetradis, and C. Wetterich, Phys. Rep. 363, 
223 (2002). 

B. Delamotte, D. Mouhanna, and M. Tissier, Phys. Rev. 
B 69, 134413 (2004). 

D. F. Litim, Nucl. Phys. B 631, 128 (2002). 
J. -P. Kownacki and D. Mouhanna, unpublished. 
M. J. Bowick et. al., J. Phys. (France) I 6, 1321 (1996). 
D. Espriu and A. Travesset, Nucl. Phys. B 468, 514 
(1996). 

Z. Zhang, H. T. Davis, and D. M. KroU, Phys. Rev. E 
48, R651 (1993). 

See the contributions of B. Duplantier in [1]. 



